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ABSTRACT 



New high resolution, high sensitivity WSRT H I synthesis observations of the spiral 
galaxy NGC 3631 are presented. In the total atomic hydrogen map, the spiral arms 
are well distinguished from the interarm regions, while the sensitivity allows detection 
of H I in all but a few isolated regions of the areas between the spiral arms. Most of 
the atomic hydrogen is located within the optical disc, but the H I extends to some 
1.5 x R op t- The H I follows the spiral arms, and streaming motions of up to ~ 15 km s _1 
(projected) can be identified from the velocity field. Assuming a constant inclination 
angle of 17°, a rotation curve is derived which is slightly falling in the outer parts of the 
disc. Analysis of a residual velocity field, obtained after subtraction of an axisymmetric 
model based on the rotation curve, confirms the existence of streaming motions near 
the spiral arms in an otherwise undisturbed disc. 

Key words: galaxies: individual (NGC 3631) — galaxies: kinematics & dynamics — 
galaxies: ISM — radio lines: galaxies — galaxies: structure 



1 INTRODUCTION 



In the study of star formation (SF) processes in discs of 
spiral galaxies, it seems natural to distinguish explicitly be- 
tween arm and interarm environments. This is however not 
usually done, partly because of the problems such an ap- 
proach entails regarding the spatial resolution of especially 
radio and molecular observations. But also in the optical 
regime, where it is easy to isolate spiral arms on images 
of nearby galaxies, photometric properties are usually stud- 
ied by means of azimuthally averaged profiles. Knapen & 
Beckman (1996) and Beckman et al. (1996) used optical im- 
ages to derive radial profiles and scale lengths for the spiral 
arm and interarm regions separately for 3 galaxies: M51, 
M100=NGC 4321 and NGC 3631, and found significant dif- 
ferences between arm and interarm profiles, with larger scale 
lengths in the arms in all three galaxies considered. For 
M100, Knapen & Beckman (1996) also included Ha, Hi, 
radio continuum and CO observations, and concluded that 
the shape of the radial (whole disc, and arm/interarm) pro- 
files is determined by SF more than by any other factor. 
Knapen & Beckman noted in particular that the Hi is en- 
hanced in the region of the star-forming spiral arms, which 
they interpret as a result, through photodissociation of part 
of the molecular gas, rather than as a cause of the SF. 

In another line of work, massive star formation efficien- 
cies (MSFE) along the spiral arms are compared directly 
with values in the adjacent interarm regions (Cepa & Beck- 



man 1990; Knapen et al. 1992, 1996). Such a comparison can 
only be made using CO, H I, and Hct data of sufficient spatial 
resolution to isolate the spiral arms. For M51 and M100, two 
galaxies where observations of these three tracers could be 
used, enhancement of the efficiencies along the spiral arms 
could be directly interpreted as evidence for triggering of 
the SF in the arms (Knapen et al. 1992, 1996). 

In order to study the role of atomic hydrogen and the in- 
terplay between gas and stars at the scale of spiral arms, H I 
observations at resolutions equivalent to at most the width 
of a spiral arm are needed. In a typical galaxy at a distance 
of ~ 15 Mpc, such as M100, M51, or in fact NGC 3631, 
this implies the need for a spatial resolution of ~ 15". Since 
the H I emission from the interarm is generally a factor 3-5 
lower than that from the arms, the need for observing the in- 
terarm H I also implies that observations of high sensitivity 
are needed. The observations described in the present pa- 
per were designed precisely to meet these goals: to reliably 
measure interarm Hi emission over the disc of NGC 3631. 

Candidate spiral galaxies for a study of the SF processes 
in their arms and disc should be objects with well-defined 
spiral arms, usually of late morphological type. They should 
also be relatively face-on objects, with say i < 30°, since pro- 
jection effects and reduced effective spatial resolution make 
it much more difficult to distinguish spiral arms in more 
inclined galaxies. This is a class of spiral galaxies not tradi- 
tionally studied in great detail through H I synthesis obser- 
vations, since the velocity information that can be deduced 
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for such low-inclination galaxies is necessarily limited. In our 
Hi study of M100 (i ~ 27°) we had to adopt and fix an in- 
clination angle when calculating a rotation curve (Knapen 
et al. 1993), and I will have to repeat that procedure in the 
present study. 

We have chosen the grand-design galaxy NGC 3631 for 
our Hi study. We follow the approach followed before for 
M100, where the H I data (Knapen et al. 1993) were used 
for detailed studies of efficiencies, SF processes, and for de- 
termining the location of the density wave resonances in the 
disc (Knapen et al. 1996; Knapen & Beckman 1996; Sem- 
pere et al. 1995, respectively). In forthcoming papers, we 
will describe similar work for NGC 3631, using the Hi data 
as presented in this paper. 

NGC 3631 was first observed at radio wavelengths by 
Roberts (1968). Other single-dish Hi observations include 
those by Fisher & Tully (1981), Tifft & Cocke (1988) and 
Staveley-Smith & Davies (1988). No synthesis observations 
have been described so far. NGC 3631 is a late-type (Sc), 
face-on spiral galaxy, looking conspicuously "normal" at op- 
tical wavelengths. It is non-barred, has no obvious compan- 
ions, and shows no other signs of important dynamical per- 
turbations. The distribution of Hn regions in the disc was 
described by Boeshaar & Hodge (1977), who also studied 
the spiral arm shape. Recently, we obtained a new high- 
quality Ha image of the galaxy, from which we catalogued 
more than 1300 individual Hn regions (Rozas, Beckman & 
Knapen 1996). 

After discussing the details of the radio and optical ob- 
servations and data reduction in Sect. 2, the distribution of 
Hi at different resolutions is described in Sect. 3. Sect. 4 is 
devoted to the kinematics, and includes a derivation of the 
rotation curve from the velocity field. The main results of 
the paper are briefly summarized in Sect. 5. 



2 OBSERVATIONS AND REDUCTION 
2.1 Atomic hydrogen 

A field centred on NGC 3631 was observed in the atomic hy- 
drogen 21 cm line with the Westerbork Synthesis Radio Tele- 
scope (WSRT) during the months of April and May, 1991. 
The total observing time of 44 hours was made up out of 
two periods of 12 h, one of 9 h, one of 8 h and one of 3 h. The 
total bandwidth of the observations was 2.5 MHz, centred 
at a heliocentric velocity of 1160 kms -1 , and divided into 
128 evenly spaced channels of 4.13 kms -1 each. The obser- 
vational parameters are fully detailed in Table 1. The data 
were reduced using the Newstar programs. Since there were 
many continuum sources within the field, 20 of them were 
subtracted from the UV data. The data were then Fourier 
transformed to a 512x512 grid of 4'.'0 x 6'.'0 (a x S) pixels. 
The resulting 128-channel datacube was Hanning-smoothed 
along the velocity coordinate, resulting in a dataset with a 
velocity resolution of 8.25 kms -1 . 

After convolving the datacube to a resolution of some 
60", 28 channels free of line emission were identified on the 
low-velocity side, and 23 channels in the high-velocity side. 
A number of channels with higher noise was discarded on 
either side. Since we are interested in structure at several 
scales and sensitivities, the data-cube was convolved from 



Table 1. Observing parameters 



Galaxy 
Telescope 

Date of observation 
Duration of observation 
Number of interferometers 
Baselines (min-max-incr) 



NGC 3631 
WSRT 

April/May 1991 

44 h 

40 

36-2700-72 



Synthesized beam (FWHM) 15'.'2 X 11'.' 2 

Pos. angle synth. beam 12? 7 

FWHP primary beam 37' 

Velocity central channel 1160.0 kms -1 (heliocentric) 

Bandwidth 2.5 MHz 

Number of channels 128 

Channel separation 4.13 kms -1 

Velocity resolution 8.25 kms -1 

Field centre (1950) ( ll h 18 m 13?2 ; 53° 26' 43") 



the original resolution of 15" 2 x 11'.' 2 to resolutions of 21" x 
14", 30" x 20", 45" x 30" and 60" x 40". The continuum was 
determined for each of these data sets by fitting a linear re- 
lation to the line-free channels, and subtracted from the line 
emission channels. The resulting five data sets were cleaned. 
Clean components were subtracted until the noise level, and 
subsequently the maps were restored by convolving the com- 
ponents with the appropriate Gaussian beam and adding the 
residuals. Five data cubes were thus produced, consisting of 
58 line channels each, at resolutions ranging from 15"2x 11'.' 2 
(hereafter full resolution) to 60" x 40". The cleaned chan- 
nel maps at 21" x 14" and 45" x 30" resolution are shown in 
Figures la and lb, respectively. Note that only every second 
channel map is shown. Noise properties of these datacubes 
and conversion factors Tt, (K)/5*(mJy) (equivalent Tj, of 1 
mjybeam -1 ) are listed in Table 2. 

The data cubes at various resolutions were used to cal- 
culate total Hi (zero-th moment), velocity (first moment) 
and velocity dispersion (second moment) maps. A careful 
inspection of the data cubes showed that such an analysis 
is in fact valid over almost the complete disc, with the pos- 
sible exception of the central region (of size approximately 
one beam), where the Hi profiles are not Gaussian-shaped 
and/or not single-peaked. I now briefly describe the pro- 
cedure followed for the full resolution data set. The first 
step was to produce a conditionally transferred data cube, in 
which values were only retained at positions where the inten- 
sity in the smoothed data cube at 30" x 20" was larger than 
2.5 times the rms noise of the smoothed maps. Pixel val- 
ues at all other positions were set to undefined. Then, noise 
peaks outside the area where Hi emission is expected were 
removed by setting pixel values at those positions to unde- 
fined. This was done interactively by inspecting the (high 
resolution) channel maps one by one, continually comparing 
with the same channel and referring to adjacent channels 
in the smoothed cube (at 30" x 20"). Finally, the resulting 
data set was used as input for the GIPSY program moments 
to calculate the total intensity, velocity and velocity disper- 
sion maps. Only emission occurring at the same position in 
at least three adjacent channels was considered true signal, 
and used for the calculation of the moment maps. 

The procedure for making the moment maps of the 
lower-resolution data sets is completely analogous. For the 
21" x 14" resolution data set, for instance, the smoothed 
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Right Aacenaion (1950) 



Figure 1. (a) Channel maps of the continuum-subtracted Hi line emission for NGC 3631, at 21"xl4" resolution. Contour levels are 
-3.8, -1.9 (dashed), 1.9 (= 2a), 3.8, 5.7, 8.6, and 11.4 mjybeam" 1 . Only every second channel is shown. The centre of NGC 3631 is 
indicated Yg' 1 MN rC ^ S 0^0 IjFP'J 011 *" vc l° c ity of each map is shown in the upper left-hand corner. The beam size is indicated as a 
natcneS ellipse in fne Dot^onwfiB corner of the last panel. 
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Right Ascension (1950) 



Figure 1. (6) Channel maps at 45"x30" resolution. Contour levels 
49.9 mjybeam -1 . Only every second channel is shown. The centre 
each map is shown in the upper left-hand corner. The beam size is 
panel. 



are -5.5, -2.8 (dashed), 2.8 (= 2<r), 5.5, 11.1, 16.6, 24.9, 33.2, and 

of NGC 3631 is indicated with a cross. The heliocentric velocity of 

indicated as a hatched ellipse in the bottom-left corner of thqTlasti 
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Table 2. Map properties 



Synthesized Beam 


rms Noise in 


Conversion factor 


(FWHM) 


channel maps 


r 6 (K)/5(mJy) 




(mjy beam -1 ) 




15'.'2 x 11"2 


0.84 


3.59 


21'.'0 x 13'.'9 


0.95 


2.09 


30'.' 1 x 19('9 


1.11 


1.02 


45'.' 1 x 29('9 


1.39 


0.45 


60" x 40" 


1.60 


0.25 



cube at 45" x 30" was used for reference. A data set of 90" x 
60" was used for reference for both the 45" x 30" and 60" x 
40" cubes, but in those cases additional restrictions were 
employed in the moment calculation, where only values of 
> 2.5a (and < — 2.5er) were used. 

2.2 Optical imaging 

Images in the B, V, R and / broad-bands were obtained in 
service time with the lm Jacobus Kapteyn Telescope (JKT) 
on La Palma, on Dec. 20, 1993. An EEV CCD chip was 
used of 1242x1152 pixels of 0"31 projected size, giving a 
field of view of around 6 arcmin. The raw images were bias- 
subtracted and fiat-fielded using dawn sky exposures, and 
photometrically calibrated using standard stars observed 
during the night. The resolution (seeing) in the reduced im- 
ages is around l'/5 . The B image as used in the present 
paper has a pixel scale of / .'62 . The position of foreground 
stars in the images was used to place the images on a cor- 
rect RA-dec grid. Since the present images were too small to 
contain a sufficient number of bright stars to warrant a good 
astrometrical solution, they were compared with the Ha im- 
age of Rozas et al. (1996), which has a larger field of view, 
and for which satisfactory astrometry could be performed 
using star positions from the HST CSC. The resulting er- 
ror in the astrometry of the optical images as used here is 
less than 0'.'2 We show the B-band image of NGC 3631 in a 
gray-scale representation in Fig. 3. 



3 DISTRIBUTION OF H I AND CONTINUUM 
EMISSION 

3.1 Global Hi properties 

The global H I line profile is shown in Fig. 2. It was produced 
by adding all the flux in each separate channel map of the 
30" x 20"resolution data set, within a square region encom- 
passing the extent of the H I emission. The systemic velocity 
(see Sect. 4) of 1155.7 kms -1 is indicated by an arrow in the 
figure. The Hi flux integral is J S dv — 51.6 ± 2 Jy kms -1 , 
in good agreement with the single-dish values given by Tifft 
& Cocke (1988) of 54.50 Jy kms^ 1 (no error indicated) , and 
Staveley-Smith & Davies (1988) of 50.7 ±4.9 Jy kms" 1 . Us- 
ing a distance to NGC 3631 of 15.4 Mpc (derived from v sys 
assuming Ho = 75kms~ 1 Mpc^ 1 ), the total atomic hydro- 
gen mass can then be evaluated as M(Hl) = 2.9(±0.1) x 
10 9 Mq. Note, however, that this value does depend on 
the true distance to the galaxy, and also on the assump- 
tion that all the atomic hydrogen is optically thin. The 
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Figure 2. Global Hi profile of NGC 3631, as derived from the 
30" X 20" resolution data set. Flux densities are corrected for pri- 
mary beam attenuation. The arrow indicates the systemic velocity 
(v ays = 1155.7 kms -1 ). 



new value compares favourably to the value reported by 
Fisher & Tully (1981) of M(Hl)= 2.9 x 10 9 M Q (corrected 
to D = 15.4 Mpc), but not to the older value of M(Hl)= 
4.5 x 1O 9 M (also at D = 15.4 Mpc) from Roberts (1968). In 
general, one can state that the integrated H I flux as derived 
from the WSRT synthesis observations agrees very well with 
previous determinations made using single-dish telescopes. 



3.2 Total H I distribution 

Figure 3 is an overlay of the 15'.'2 x 11'.' 2 total Hi (zero-th 
moment), or column density of hydrogen, map of NGC 3631 
on the B-band CCD image. The H I generally traces the 
spiral arms, although the correspondence between optical 
and Hi features is not unique. The Hi extends further out 
than the optical disc, especially in the SW region of the H I 
extension (see below and Fig. 4b), but it does not do so in 
the dramatic fashion as seen in e.g. M101 (Allen et al. 1974; 
van der Hulst & Sancisi 1988) or NGC 628 (Kamphuis & 
Briggs 1992), two galaxies of morphological type comparable 
to NGC 3631. 

Figures 4a and 4b show the Hi column density distri- 
bution (zero-th moment maps) in NGC 3631 at resolutions 
of 15" 2 x 11" 2 and 45" x 30", respectively. The higher res- 
olution map (Fig. 4a) shows that the Hi is mostly concen- 
trated in the spiral arms, but that Hi is detected all over 
the disc, including the interarm zones and the central region 
of the galaxy. At 45" x 30"resolution (Fig. 4b) the general 
arm pattern can still be recognized but is mostly lost in the 
lower resolution. An interesting feature is the extension on 
the SW side of the disc, at low levels. This is reminiscent 
of the distribution of Hi in M100, but is not accompanied 
by a kinematic signature, as is the case in M100. Similar H I 
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Figure 3. 15" 2 X 11" 2 total Hi column density map (contours) overlaid on a grayscale representation of a S-band CCD image of 

NGC 3631. Contour levels are 2.0, 4.1, 6.1, 8,2, 12.2 and 16.3 X 10 20 atoms cm . Beam size is indicated. Note that the local minima in 

the Hi distribution can be recognised by comparison with the grey-scale representation of the same Hi map in Fig. 4a. n n 
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Figure 4. a (upper) Total Hi (column density) map of NGC 3631 at 15" 2 X 11'' 2 resolution. Contour levels as in Fig. 3. The centre is 
indicated with a cross, and the beam size is indicated in the lower left corner. Outermost contour is interrupted by undefined points, b 
(lower) As Fig. 4a, now at 45" X 30" resolution. Contour levels are from 0.85 to 10.37 in steps of 1. 19 xlO 20 atoms cm . Beam size is 
indicated. 
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o ' 

12 3 4 

Radius (arcmin) 

Figure 5. Radial Hi surface density profile of NGC 3631, as 
derived from the 30" X 20" total Hi map by averaging in elliptical 
rings of fixed inclination and position angles. The extent of the 
optical disc (R op t) is indicated. Note that the Hi outside R ~ 3' 
resides in the SW extension, rather than being distributed in the 
disc (see text). 

asymmetries were first recognised by Baldwin, Lynden-Bell 
and Sancisi (1980), and seem to be a frequent phenomenon. 

3.3 Radial H I profile 

A radial H I surface density profile was derived from the 
30" x 20" total H I map by averaging in elliptical rings with 
fixed inclination and position angles (i — 17° and PA= 150° 
respectively, see Sect. 4), and is shown in Fig. 5. The extent 
of the optical disc (R op t = (0.5 x D25); from de Vaucoulcurs 
et al. 1991) is indicated in the Figure. The profile shows 
clearly some features that were already obvious from the 
overlay of the total Hi map on the B-band optical image 
(Fig. 3), namely that the Hi disc is not a whole lot more 
extended than the optical disc, and that most of the Hi sits 
well inside the optical disc. The radial Hi profile shows a 
central depression, not uncommon at all for both barred and 
non-barred galaxies (see e.g. Broeils & van Woerden 1994). 
The profile peaks at a radius of just over 1 arcmin, and falls 
off rapidly after that radius, coinciding with the end of the 
region of the star forming spiral arms as seen in the optical 
image. This situation is reminiscent of that in M100, where 
also the Hi is enhanced in the region of the SF spiral arms. 
Knapen & Beckman (1996) interpret this enhancement in 
M100 as a result of the SF activity, leading to photodisso- 
ciation of part of the molecular gas and to the production 
of Hi, rather than as the origin of the SF. Careful compari- 
son with especially molecular gas observations is needed to 
confirm that a scenario of Hi production as a result of SF 
activity is in fact the preferred one also for NGC 3631. 

The radial Hi profile extends out to around R — 4', 
but inspection of especially Fig. 4b shows that the H I only 
extends to this radius in the SW, elsewhere the disc can not 
be defined for radii larger than R ~ 3'. This implies that 
the surface density beyond 3 arcmin in Fig. 5 is no longer a 
ring average, as in the inner parts, but is due mostly or even 
exclusively to the H I in the SW extension. Comparison of 




RA 



Figure 6. 15'.'2 X 11'.' 2 21 cm continuum map (contours and 
grayscales). Contour levels are —1.0, —0.5 (dashed), 0.5 (= 2a), 
1.0, 1.5, 2.0, 2.5, 3.3 and 4.0 mjybeam- 1 . Beam size is indicated. 

the whole radial H I profile of NGC 3631 with those for other 
Sc galaxies (Broeils & Rhee 1996) shows that the H I disc 
in NGC 3631 is rather small. For NGC 3631 we find a value 
for the ratio of Hi to B diameter DH1/D25 of 1.1 (using 
am = 1 Mq pc -2 for the definition of Dhi)- Broeils & Rhee 
(1996) find that these ratios for the around 15 Sc galaxies in 
their sample are between 1.1 and 2.5, placing NGC 3631 at 
the lower extreme. The peak surface density in NGC 3631, of 
some 7Mq pc -2 , and the average surface density of around 
4.5 Mq pc -2 , are quite normal when compared to both the 
Sc galaxies, and also to the galaxies of other morphological 
types, from Broeils & Rhee (1996). 

3.4 21 cm Continuum 

The 21 cm continuum map, at 15"2 X 11'.' 2 resolution, is 
shown in Fig. 6. This is the image that was subtracted from 
the data cube of the same resolution. As described in Sect. 2, 
a number of continuum sources were deleted from the data 
set at the stage preceding the making of the channel maps, 
but still a number of grating rings are present in the final 
continuum map. This makes interpretation of the features in 
the map, and especially the determination of fluxes, very dif- 
ficult. What can be deduced from the map is that the 21 cm 
continuum emission from NGC 3631 is centrally peaked, and 
follows the spiral arm shape. No obvious radio continuum 
point sources are seen in the disc, thus no recent supernova 
events can be identified. 



4 KINEMATICS 

The velocity field (first moment map) of NGC 3631 at 21" x 
14" resolution is shown in Fig. 7a in a contour plus grayscale 
representation. Similar maps were in fact produced at all 
resolutions; only one is shown in the present paper. The 
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Figure 7. a. (upper left) Hi velocity field of NGC 3631 at 21" X 14" resolution. Contour levels are from 1100 kms -1 to 1200 kms -1 
in steps of 10 kms , where the lower values are found to the SE, and the first white contour is at 1160 kms" 1 . Grayscales indicate 
rourfdg th^samc^mr^^ ^n vchliinttes. Beam size is indicated, b. (upper right) Model velocity field as determined from the rotation curve 
^cC^cxt^^Con^our^ ana gray l(«as as in Fig. 7a. c. (lower left) Residual velocity map, obtained by subtracting the model (Fig. 7b) from 
the velocity field (Fig. 7a). Contours are at —10, —6 and —2 kms - (black) and 2, 6, and 10 kms - (white), with grayscales indicating 
the same range and higher values coded darker, d. (lower right) Position- velocity diagram along the major axis (<f> = 150°) of the 21" X 
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velocity field in Fig. 7a shows a generally regular shape, but 
a few features are worth pointing out explicitly. The shape 
of the isovelocity contours indicates that the position angle 
of the major axis is practically constant over the disc. The 
closed contours representing more extreme velocity values 
point out a constant or slightly falling rotation curve in the 
outer half of the disc, shown to be indeed the case in the 
next paragraph. At several positions in the disc, notably 
about 1' W and ~ 2' E of the nucleus, and showing up most 
clearly along the minor axis, are deviations from the regular 
shape of the isovelocity contours that can be recognized as 
streaming motions due to a density wave near the H I spiral 
arms (e.g. Rots 1975). From the displacement of the contours 
one can estimate an amplitude of ~ 15 kms -1 along the 
line of sight, or ~ 50 kms -1 in the plane of the disc, after 
deprojection, assuming an inclination angle of 17° (RC3, 
see also next paragraph). Such values are somewhat large, 
but not outside the range of values found in other galaxies 
for density wave streaming motions (e.g. Visser [1980] for 
M81, Rots et al. [1990] for M51, and Knapen et al. [1993] 
for M100). 
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4.1 Rotation curve 

I derived the rotation curve from the velocity field (first 
moment maps) at different resolutions, using a procedure 
described in detail by Begeman (1989). In this procedure, 
the galaxy is divided into a set of concentric rings, each of 
which is described by a set of parameters i (inclination an- 
gle), <f> (position angle of the major axis), and v c (rotational 
velocity) . Additional parameters are the centre of each ring, 
and the systemic velocity v sys of the galaxy. These param- 
eters are fitted using a least squares algorithm. Data points 
within each ring are weighted by | cos(t9)|, where 8 is the az- 
imuthal angle from the major axis. Data points within 20° 
from the minor axis are excluded from the fits. 

NGC 3631 is a galaxy with a low inclination angle 
(i — 17°, de Vaucouleurs et al. 1991) which makes it im- 
possible in practice to fit both v c and i. Rather than as- 
suming a constant rotational velocity throughout the disc, 
a constant inclination angle was used. Using the value of 
i — 17° (RC3) results in values for the rotational velocity of 
v c ~ 140 kms -1 , which agree quite well with the synthetic 
rotation curve given by Rubin et al. (1985) for a galaxy like 
NGC 3631. Rubin et al. give a value of v ro t = 132 kms -1 
for an Sc galaxy of Mb = —20 mag, the absolute magnitude 
of NGC 3631 assuming D = 15.4 Mpc and the value for mj 
from the RC3. The agreement is surprisingly good given the 
high uncertainty in the determination of an inclination angle 
of only 17° . Note that the use of a different value for i does 
not change the shape of the rotation curve, but only scales 
the values for v c (larger v c for smaller i). 

Starting out with the 60" x 40" velocity field, I first 
fitted the position of the dynamical centre as q(1950) = 
n h lg m 12=75^0:07), 5(1950) = 53° 26' 43'.'0(±l") by fixing 
i, <f> and v sys at reasonable values. This position was deter- 
mined from the part of the disc at radii larger than 40", 
where the errors in the fit were acceptable. The dynami- 
cal centre coincides with the optical centre of the galaxy (as 
given in the RC3) in declination, but is offset from it by 0.55 
seconds of time, or some 5 arcsec, in right ascension. As a 
second step, the position of the centre and the inclination 



Figure 8. Rotation curve for NGC 3631 (upper panel): whole 
disc (full dots), receding side only (open circles), and approach- 
ing side (triangles). Formal errors from the least-squares fit arc 
indicated, but the uncertainty due to the low inclination of the 
galaxy is larger. Fitted position angle of the major axis (N over 
E) is shown as a function of radius in the lower panel. For both 
panels, points from the centre out to 2'.4 are derived from the 
21" X 14" resolution velocity field, points at R > 2.'5 are from the 
60" X 40" velocity field. 



angle were fixed, which resulted in a satisfactory fit to the 
systemic velocity of the galaxy, v sys = 1155.7 ± 0.8 kms -1 . 
This is in good agreement with the values determined from 
single-dish Hi observations (e.g. v sys — 1156 ± 2 kms -1 , 
Staveley-Smith & Davies 1988). 

The final rotation curve fit for the whole disc was now 
made fixing the position of the centre and v ays at the val- 
ues described above, i at 17°, and fitting cj> and v c at ra- 
dial points with 20" (similar to half beam) spacing. Fits for 
the receding and approaching halves of the galaxy were also 
made. 

Rotation curves for the whole disc, and for receding and 
approaching halfs separately, were then fitted for all other 
resolutions, using as width of the rings always half the (mi- 
nor axis) beam size. Following exactly the same procedure as 
described for the 60" X 40" velocity field, values for the dy- 
namical centre and systemic velocity were first checked, and 
<j} and v c subsequently fitted for all the other cubes at lower 
resolution. The resulting final rotation curve for NGC 3631 
is shown in Fig. 8. This Figure is a composite of the fits 
to the 21" x 14" resolution velocity field (for radii out to 
2!4 ) and to the 60" x 40" velocity field (larger radii). The 
top panel shows the rotation curve for the whole disc, and 
approaching and receding sides separately, while the lower 
panel shows the run of the fitted position angle of the major 
axis (measured N over E) against radius. It is clear that <f> 
is essentially constant over most of the disc, at a value of 
4> = 150° ±2°. It reaches significantly larger values (up to 
4> = 160°) only around R = 55", which is probably caused 
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by the streaming motions prevalent there (see previous sec- 
tion). 

The rotation curve rises slowly out to R ~ 1.'4, after 
which it declines somewhat. This decline does however de- 
pend critically on the constancy of the assumed inclination 
(a chance of only 1 or 2° would suffice to cause the ob- 
served drop). Over most of the disc the rotation curves for 
approaching and receding sides of the disc behave very much 
like the total rotation curve. Apart from the outer points, 
the only deviations are found in the inner Of 6. These may 
well be caused by the adopted central position, which was 
determined mostly from the outer regions of the disc, and 
from data at 60" x 40" resolution. An overlay of the (21" x 
14") rotation curve on a position- velocity diagram along the 
major axis (Fig. 7d) confirms that the derived rotation curve 
is very reasonable, even though the inclination of the galaxy 
is low. 

4.2 Model velocity field 

The fitted runs of v c and (f> for the 21" x 14" data (whole disc 
fit; see previous section) were used to construct an axisym- 
metric model velocity field. All other parameters (v sy s,i and 
central position) were kept at the values described above. 
The resulting model velocity field is shown in Fig. 7b. The 
counterclockwise deviation of the isovelocity contours near 
the outer edge of the model is due to a pair of high val- 
ues of (f>, which may be artifacts. The model looks gener- 
ally smooth, although it is clear that the streaming mo- 
tions discussed before show up to a certain degree. The 
model was subtracted from the velocity field to produce 
a residual velocity map, which is shown in Fig. 7c. Values 
for the residual velocity of \v rE s\ > 15 kms -1 occur exclu- 
sively in the outer 40" of the map. Near the centre, where 
large velocity gradients occur, residual velocities are posi- 
tive, around v res — 10 km s _1 . In the disc, the residual veloc- 
ity field shows a continuous region of negative values (dark 
in Fig. 7c) , aligned more or less along the direction of the mi- 
nor axis. This feature may well result from the fact that data 
points within 20° from the minor axis were excluded from 
the fit leading to the rotation curve. The residual velocity 
map also shows the spiral arms, in the form of positive resid- 
ual velocities where streaming motions can be identified in 
the velocity field (Fig. 7a). This is most obvious in the form 
of the curved feature in Fig. 7c starting N of the nucleus, 
curving to the left (E) and then down (S), although other 
features in the map can also be identified to lie near or on 
the spiral arms (compare Fig. 3). The absence of a symmet- 
ric pattern surrounding the centre confirms that this galaxy 
does not have a dynamically important non-axisymmetric 
component (such as a bar). 



5 SUMMARY 

In this paper, I describe new WSRT H I aperture synthe- 
sis observations of the grand-design late-type spiral galaxy 
NGC 3631. The highest spatial resolution of the produced 
data set is of 15'.'2 x 11" 2, which allows resolving the spiral 
arms from the interarm disc, but data cubes at several lower 
resolutions were also produced, with the lowest resolution 
described here being 60" x 40". These Hi data will be used 



in subsequent work for detailed studies of SF processes in 
and outside the spiral arms, and of the spiral structure. The 
main parameters for NGC 3631, as taken from the literature 
or determined in the present paper, are listed in Table 3, and 
the main results can be summarized as follows: 

(i) H I is detected all over the disc of the galaxy, with good 
detections, even at the highest resolution, in the interarm 
regions and in the centre. The H I generally follows the spiral 
arm pattern as seen in optical images. 

(ii) The Hi extends to about 1.5 x R op t, but most of the 
H I is found well within the optical disc. The radial H I profile 
peaks just outside R = 1', in the region of the star-forming 
spiral arms. The profile shows a slight central depression. 

(iii) Streaming motions due to density waves can be iden- 
tified near the locations of the spiral arms in the high- 
resolution velocity field. They have maximum amplitudes of 
~ 15 kms -1 , or ~ 50 kms -1 when assuming an inclination 
angle of 17°. No significant other deviations from axisym- 
metry are seen in the velocity field. 

(iv) A rotation curve as derived from the velocity field 
rises slowly to its maximum at R ~ l!4, after which it de- 
clines somewhat. The exact values of the rotational velocity 
v c could not be determined due to the ambiguity between 
the angle of inclination and v c , but assuming i — 17° we find 
that v c ~ 140 kms -1 , in general agreement with synthetic 
rotation curves for a galaxy like NGC 3631 (Rubin et al. 
1985). The position angle of the major axis <j> — 150° ± 2° 
is practically constant over the disc. 

(v) After subtracting a model velocity field, obtained 
from the rotation curve, from the original velocity field, the 
residual map confirms the existence of streaming motions in 
the H I, and the absence of significant other disturbances in 
the disc. 
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Table 3. Parameters and results for NGC 3631 



Parameter 


Value 


Notes 


Morphological type 


SA(s)c 


1 


Optical centre (1950) 


ll h 18 m 13?3 , 53° 26' 43" 


1 


Dynamical centre (1950) 


ll h 18 m 12!75 (±0?07), 53° 26' 43'.'0 (±1") 


2 


Systemic velocity (heliocentric) 


1155.7±0.8 kms -1 


2 


Distance 


15.4 Mpc 


3 


m B 


10.97 ±0.14 


1 


Optical size (0.5 X D25) 


2.'5 


1 


H 1 radius 


4.'0 ± 0^2 


2 


Inclination 


17° 


1, 4 


Position angle of major axis 


150°±2° 


2 


H 1 flux integral 


51.6 Jy kms" 1 


2 


Total atomic hydrogen mass 


2.9 x 10 9 M Q 


2 



note 1: RC3 (de Vaucouleurs et al. 1991) 
note 2: This paper 

note 3: From the systemic velocity, assuming Hq = 75 kms" 1 Mpc" 
note 4: Adopted (see Sect. 4) 
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